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ABSTRACT 

Gamma ray bursts have been divided into two classes, long-soft gamma ray burst and short-hard 
gamma ray burst according to the bimodal distribution in duration time. Due to the harder spectrum 
and the lack of afterglows of short-hard bursts in optical and radio observations, different progenitors 
for short-hard bursts and long-soft bursts have been suggested. Base d on the X-ray afte rglow obser- 
vation and the cumulative red-shift distribution of short-hard bursts, iNakar et al.l (|2006[ ) found that 
the progenitors of short-hard bursts are consistent with old populations, such as mergers of binary 
neutron stars. Recently, the existence of two subclasses in long-soft bursts has been suggested after 
considering multiple characteristics of gamma-ray bursts, including fluences and the duration time. 

In this work, we extended the analysis of cumulative red-shift distribution to two possible subclasses 
in L-GRBs. We found that two possible subclass GRBs show different red-shift distributions, especially 
for red-shifts z > 1. Our results indicate that the accumulative red-shift distribution can be used as 
a tool to constrain the progenitor characteristics of possible subclasses in L-GRBs. 
Subject headings: gamma rays: bursts — binaries: close — supernovae: general 



1. INTRODUCTION 

Traditionally, gamma-ray bursts (denoted as GRBs) 
have been divided into two classes which are separated 
by the duration time of ga mma-ray bursts Tgg = 2 
sec ()Kouveliotou et all 1X9931 ). Since the short-duration 
gamma-ray bursts had relatively hard observed spec- 
tra than long-duration gamma-ray bursts, two classes 
are usually called as long-soft gamma ray burst (de- 
noted as L-GRB) and short-hard gamma-ray burst (de- 
noted as SHB). From the afterglow observations of L- 
GRBs, the association between the L-GRBs and su- 
pernovae/hypernovae has been suggested. These ob- 
servations suggest that the progenitors of L-GRBs are 
massive giant stars which are collapsin g. One of the 
best candidate is the Collapsa r model (jWooslevI 1X9931 : 
iMacFadven and WooslevI 1X9990 in which rapidly rotat- 
ing black hole with Kerr parameter about 0.8 is formed 
after GRB. Recently, spinning black holes with high Kerr 
parameters, a^, « .8, have been observed in soft X-ray 
black hole binaries (jMcChntock et al.ll2006HShafee et all 
12006'). These observations support the idea that the 
rapidly rotating black holes in soft X-ray black hole bi- 
naries are the r emnants of L-GRBs ()Brown et al.ll20dol : 
iLee et al]|2002[ ). providing the natural sources for Col- 
lapsars. 

Due to the differences in the observational character- 
istics of SHBs compared to those of L-GRBs, such as 
lack of afterglow observations in optical and radio band 
and harder spectrum, etc., different progenitors for SHBs 
have been suggested. Mergers of binar y neutron stars 
(|Eichler et al.l 1X9891 : iNaravan et~aLllX992D are among the 
best candidates. After recent X-ray afterglow observa- 
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tions by Swift and HETE-II, the associations between 
SHBs and the host galaxies were reported (see Table [T] 
and references there in). Based on the X-ray afterglow 
observation of SHB s and the cumulat ive red-shift dis- 
tribution of SHBs, INakar et all ()2006l ) found that the 
progenitors of SHBs are consistent with old populations. 
This finding supports mergers of binary neutron stars as 
progenitors of SHBs. 

Recently, the existence of two subclasses in L-GRBs 
has been suggested after considering both fluen ces and 
the duration time (jGhattopadhvav et al.ll2007[) . They 
divided L-GRBs into Clusters II and HI which were sep- 
arated by the fluences combined with the duration time. 
They classifled L-GRBs with higher fluences as Cluster 
HI and suggested that their progenitor s are the mas- 
sive collapsing stars, such as C ollapsar (jWooslevI IX993l : 
IMacFadven and WooslevI [19990 . As a possible origin of 
Cluster II GRBs, they suggested neutron-star, white- 
d warf binaries. If there exist t wo subclasses in L-GRBs 
as lChattopadhvav et al.l ()2007| ) suggested and only Clus- 
ter III GRBs are associated with massive collapsing gi- 
ants which are the progenitors of supernovae or hyper- 
novae, the observed red-shift distribution of two sub- 
classes might show different behaviour. 

In this work, we extended the work of INakar et al.l 
(l2006l ) to two possible subclasses in L-GRBs. We found 
that Clustters II and HI show different red-shift distribu- 
tion at high red-shift z > 1. Our results indicate that the 
accumulative red-shift distribution can be used as a tool 
to constrain the lifetimes of possible subclasses in GRBs. 
In Sec. 121 the numerical methods which we used in our 
analysis are summarized (jNakar et al.ll2006l ). The esti- 
mated lifetime of SHBs are summarized in Sec. |3l We 
confirmed that the progenitors of SHBs arc consistent 
with old population with lifetime t,^ = 6.5 Gyr, indicat- 
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ing that SHBs are from the mergers of binary neutron 
stars. In Sec. IH the cumulative red-shift distributions 
of two possible subclasses in L-GRB have been summa- 
rized. Our results show that there are clear differences 
in the red-shift distributions of two possible subclasses 
in L-GRBs for the red-shift z > 1. This suggests that 
the cumulative red-shift distribution can be used as a 
tool to distinguish the subclasses in L-GRBs. Our final 
conclusion follows in Sec. [5l 

2. GAMMA-RAY BURSTS OBSERVATION RATE 

To examine the progenitor's life time and the local 
rat e of GRBs, we us ed the local GRB rate suggested 
bv lNakar etatl ()2006f) . 



SFR(z/)/(t(z) 



/ •,^ dt , 



(1) 



where -Rgrb(z) is the intrinsic GRB rate (per unit co- 
moving volume and comoving time), SFR(z) is the local 
star formation rate (SFR) at red-shift z (per unit co- 
moving volume and comoving time), t{z) is the age of 
universe at red-shift z, and /(r) is the fraction of GRB 
progenitors that are born with a lifetime r. Note that 
the absolute value of i?GRB is not important because we 
normalized the final cumulative distribution to unity. In 
this work, we used three different parameterizations of 
the global st ar formation rate per un it c omoving vol- 
ume foUowing lPorciani &: Madaul (|2001h and lNakar et al.l 
((2006.) . 
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where ii{z^ in the s tandard ACDM Universe model 
(jCopeland et al.ll2006[ ) is 
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Fig. 1. — Normalized star formation rate (SFR) as a function 
of red-shift z. Due to the uncertainty in the estimation of SFR at 
high red-shift, z > 1.5, th ree different choices of par ametrization, 
Eq. are used followine lPorciani fc Madaul 11200 n't . 



with Hubble constant Hq. In Fig.[l] we summarized the 
normalized SFR. One can see the clear differences o f SFR 
at high red-shift z > 1.5 ([Porciani fc Madaul lMll). 

We used two diff erent forms of /(r) in Eq. ^ following 
iNakar et al.) (|2006l ). 
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The first is the power-law distribution with which the 
dominant portion of GRBs were produced at the early 
evolutionary stage of their progenitors. The latter is 
the log-normal distribution in which r* is the mean time 
and <j is the di spersion. In this work we used a = 0.3 
(|Nakar et all 



In order to compare the local GRB birth rate at red- 
shift z with the current observation on earth, we used 
the differential GRB observation rate dNobs/dz] 
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where the luminosity distance dL{z) is 
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dL{z) = (l + z) 
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and the k{z) correction, which inc ludes the conversi on 
from energy flux to photon flux, is ()Nakar et al.ll2006l ) 

k{z) « 2 X 10"'^(1 + zy^-"^ ergs. (8) 

S{P) is the probability for detecting the bursts with a 
peak photon flux P which in turn depends on the lumi- 
nosity L and the red-shift z as well as on the spectrum 
bursts. Since the exact form of S{P) is unknown, in this 
work we used the simplest form S{P) — 1. The parame- 
ter (3 in Eq. ^ is the power in the luminosity function 
(/)(L) = (f>oL~^ and we used = 2 wh ich is the best fit 
parameter found bv lNakar et al.l (|2006D . 

3. LIFETIME OF SHORT-HARD BURSTS 

The observed red-shifts of SHBs are summarized in 
Table [T] For GRB 050813, two values of red-shifts are 
given due to the uncertainty in the ass ociation between 
the putative host galaxy and the SHB ()Prochaska et al.l 
[2005; Bcrger 2006). In our analysis, we used 7 SHBs 
without GRB060502B in which the association with the 
putative host galaxy is uncertain. Our final conclusion on 
the lifetime of SHBs remains the same even if we include 
GRB060502B. 

The cumulative red-shift distribution of SHBs are sum- 
marized in Fig. [21 From the upper panel in this figure, 
one can see that it's very hard to fit the observed distri- 
bution using power-law distribution. We also confirmed 
that SFR2 and SFR3 produce similar results. From these 
results, one can confirm that SHBs are not associated 
with explosive phenomena of giants, such as supernovae 
or hypernovae. In the lower panel of Fig. [21 the log- 
normal distribution Eq. ([5]) with various mean time 
is considered. We found that the old population with 
~ 6.5 Gyr gives the best fit to the observation. This 
results indicates that SHBs are associated w ith the merg- 
ers of neutron star or black hole binaries (|Nakail [20071 : 
iBethe et alll2007HLee et al.ll2007[ ). 
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SHB 


Red-shift (z) 


Type 


Association^ 


Ref. 


050509B 


0.225 


E(c) 


3 — 4cr 


1 


050709 


0.16 


Sb/c 


Secure 


2 


050724 


0.258 


E/SO 


Secure 


3 


050813 


0.72 or 1.8 


E/SO 




4 


051221A 


0.546 


Late 


Secure 


5 


060502B(?) 


0.287(?) 


Early 




6 


790613 


0.09 


E/SO 


So- 


7 


000607 


0.14 


Sb 


2u 


7 



TABLE 1 

Putative host galaxy properties of S HBs (INakarII200 
References: IVGehrels et al.' ("2005'): 'Fox et al.' f200 
[Bloom et .'\l 



(2006); 2) Fox et al. (2005); Prochaska et al.I 
()200 5); CoviNO et al. (2006); 3) Berger et al. (200^; 4) 
Prochaska et al. (2005); Berger (2006); 5) 
[Soderberg et al. (2006); 6) Bloom et al.. (i2007l) ; 7) 
I Gal- Yam et a l. (2005). f The confidence level of the 
association between the putative host and the SHB. 
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Fig. 2. — The cumulative red-shift distribution of SHBs for two 
different choices of /(r). Here we used SFRl and f2m = 0.3 for the 
matter contents of the Universe. The upper panel corresponds to 
power-law distribution /(r) = and the lower panel corresponds 
to the log- normal distribution, Eq. ((Sj. Thick solid line corresponds 
to the observed data. Two different curves are the consequence of 
the uncertainty in the red-shift estimation of SHB050813. 



Fig. 3. — The cumulative red-shift distribution of SHBs with 
SFRl for different choices of Urn in ACDM model. Two different 
mean values of delay time, 6.5 Gyr and 3 Gyr, are used for the 
comparison. 

In previous results, we used the matter content in the 
Universe il„i — 0.3 which is close to the current observa- 
tion. We investigated SHB progenitor life time by chang- 
ing the matter contents Qm in Eq. ([3]). The results are 
summarized in Fig. [3] in which we considered two differ- 
ent choices of the mean time r*. From Fig. [3l one can 
see that the old population with mean time r* = 6.5 Gyr 
is consistent with the observation independently of the 
matter contents. We also tested other Universe models, 
such as Phantom CDM model, and found that our con- 
clusion is independent of the particular Universe models. 
This observation is consistent with the fact that most of 
the observed SHBs are distributed in low red-shifts. 

4. SUBCLASSES IN LONG-SOFT GAMMA-RAY BURSTS 

Due to the earlier afterglow observations, the progen- 
itors of L-GRBs are believed to be related with those 
of supernovae or hypernovae. However, recently the 
pos sibility of two sub-cla s ses in L-GRBs was suggested 
by IChattopadhvav et all (|2007l ). They divided long- 
duration GRBs into Cluster 11 for low fluence and Cluster 
HI for high fluence, which were separated by the line, 

1.6 X 10~V790 ergs/cm^ (9) 
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Fig. 4. — The cumulative red-shift distribution of L-GRBs. Thick 
Unes correspond to the observed distributions. The thiclc soUd Hne 
corresponds to the 43 Cluster III GRBs and the thick dotted line 
corresponds to the 48 Cluster II GRBs. Two clusters show clear 
difference in the red-shift region 2 > 1. Thin lines correspond to 
the lowest cumulative distribution for given star formation rates 
with /3(r) = 

where Ft is the fluence and Tgo is the time corresponding 
to 90% of the flux. They showed that Cluster II GRBs 
have nearly constant isotropic energy output of 10^^ ergs 
while Cluster III GRBs have isotropic energy in wide 
ranges from 10^^ to 10^'' ergs. As possible origins, they 
suggested mergers between neutron star and white dwarf 
for Cluster II and collapse of massive stars for Cluster III. 

In this work, we investigated t he cumulative red - 
shift distribution of Cluster III using iNakar et al.l (pOOl ) 
method. If Cluster III is related with the collapse of mas- 
sive stars, their cumulative distribution should be able to 
be described by Eq. ([4]). The results are summarized in 
Fig. m In this figure, we tested two different samples of 
L-GRBs with known red-shifts, 48 Cluster II (thick dot- 
ted line) and 43 Cluster III (thick solid line) which were 
divided by Eq. Three thin lines are the lowest cu- 
mulative distributions for given star formation rates. In 
Fig.m one can see that independently of SFR models the 



Cluster II L-GRBs cannot be described by the power-law 
distribution for the red-shift z > 1. On the other hand, 
even though it's too early to make any firm conclusion, 
distributions with only Cluster III can be easily described 
by the power-law distribution. This might indicate that 
there exist two subclasses in the long-duration GRBs. 
Main deviations between two clusters in Fig. 2] are in the 
red-shift region z > 1. 

5. CONCLUSIONS 

In this work, by exteiiding the work of Nakar et al. 
(|Nakar et al.ll2006l : [N"akaill2007f) . we investigated the con- 
straints on the progenitor lifetime of SHBs. We con- 
firmed that the SHBs are consistent with old population 
with mean time r, = 6.5 Gyr. We also confirmed that 
this conclusion is independent of the details of the cos- 
mology models. This results support the conclusion that 
the progenitors of SHBs are consistent with old popu- 
lation, such as merging compact star binaries (neutron 
star-ne utron star binaries or neutron star-black hole bi- 
naries ) (INakar et al.ll2006l : lNakad[2007t iBethe et al.ll2007l : 
iLee et al.ll2007E 

We also investigated two subclasses Cluster II and 
HI L-GRB s ( Chatto padhvav ct al. 2007) using the same 
method of iNakar et all ([20061). We found that the power- 
law distribution is more consistent with Cluster HI GRBs 
than with Cluster II GRBs. We also found that the cu- 
mulative red-shift distribution of two clusters show clear 
difference in the red-shift region z > 1. This result sup- 
ports the existence of two subclasses in L-GRBs. How- 
ever, in order to draw firm conclusions, many effects 
which are not included in our analysis, e.g., the con- 
nection between GRBs and the metallicity of the host 
galaxy, different characteristics of cosmological GRBs 
and subluminous GRBs, observability premium for dif- 
ferent type of GRB progenitors, etc., have to be consid- 
ered. 
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